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Abstract. We consider the Cauchy problem for (energy-subcritical) nonlin- 
ear Schrodinger equations with sub-quadratic external potentials and an addi- 
tional angular momentum rotation term. This equation is a well-known model 
for superfluid quantum gases in rotating traps. We prove global existence (in 
the energy space) for defocusing nonlinearities without any restriction on the 
rotation frequency, generalizing earlier results given in [11, 12]. Moreover, we 
find that the rotation term has a considerable influence in proving finite time 
blow-up in the focusing case. 



1. Introduction 

Ever since the realization of Bose-Einstein condensation (BEC) in dilute atomic 
gases, much attention has been given to dynamical phenomena associated to its 
superfluid nature. One remarkable feature of a superfluid is the appearance of 
quantized vortices, cf. [1] for a broad introduction to these kind of phenomena. In 
physical experiments, the BEC is thereby set into rotation by a stirring potential, 
which is usually induced by a laser [19, 20, 21, 24] (see also [4] for numerical simula- 
tions). The corresponding mathematical model is a nonlinear Schrodinger equation 
(NLS) with angular momentum rotation term, i.e. 

(i.i) ^ t v> = -yAV> + A|^|^ + F(x)v>-ft-£V>, (i,i)eixi 3 , 

where ip = tp(t,x) is the complex- valued wave function of the condensate. In the 
physics literature, (1.1) is known as the Gross- Pitaevskii equation for rotating Bose 
gases. The coupling constant A £ 1 can be experimentally tuned to account for 
both defocusing (A > 0) and focusing (A < 0) nonlinearities. The potential V(x) 
describes the magnetic trap and is usually assumed to be of the form 
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(1.2) V(x) i^-- 
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Finally, il ■ L denotes rotation term, where 
(1.3) L:=-ixAV 

is the quantum mechanical angular momentum operator and 6 R 3 is a given 
angular velocity vector. For a rigorous derivation of (1.1) in the stationary case, 
we refer to [15]. Furthermore, we remark that the appearance of quantized vortices 
has been rigorously proved in [23] , by means of a spontaneous symmetry breaking 
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for the ground state of the stationary equations (provided A > is sufficiently big). 
For further mathematical results in this direction we refer to [1] and the references 
given therein. 

The aforementioned works illustrate the fact that there is a considerable amount 
of mathematical studies devoted the stationary equation. On the other hand, the 
time-dependent equation (1.1), has not been given as much attention, even though 
it is considered to provide the basis for a dynamical description of vortex creation. 
Indeed, except for numerical simulations [4] , we are only aware of [11, 12] providing 
rigorous results for (1.1). In [12] global well-posedncss of the Cauchy problem (in 

the energy space) is proved in the case where A > 0, V(x) = -7f\x\ 2 7 i.e. an isotropic 
confinement, and |0| = 7. The analogous result for in d — 2 spatial dimensions is 
given in [11]. 

Remark 1.1. Let us mention that in [16, 17], the NLS model (1.1) is also rigorously 
studied. The results, however, mainly concern an asymptotic regime, the so-called 
semi-classical limit, and are thus very different from the present work. 

In view of these results the main goal of our work is twofold: First, we shall prove 
global well-posedness of (1.1) in the defocusing case, without any restriction on \Q\ 
or The latter is needed to describe actual physical experiments, which 

often require \Q\ ^ 7. To this end, we shall show that by a suitable time-dependent 
change of coordinates, we can transform equation (1.1) into a nonlinear Schrodingcr 
equation without rotation term but with a time- dependent trapping potential. In 
a second step, we shall analyze the possibility of finite time blow-up of solutions, 
in the case of a focusing nonlinearity. Recall that finite time blow-up means, that 
there is a T* < +00, such that 

ffin ||VV(*)|k2 =+00. 

As we shall see, the usual proof of finite time blow-up, based on the classical virial 
argument of Glasscy (sec e.g. [8]), in general does not go through in a straight- 
forward way, due to the influence of the rotation term. Instead, it has to be slightly 
modified, yielding blow-up conditions which depend on |Q|, and which coincide with 
the usual conditions in the limit |f2| — > 0. 

2. Mathematical Setting and Main Result 

In the following we shall consider the Cauchy problem for the following, slightly 
more general NLS type model 

(2.1) idt<l>= -\A^ + \\tf;\ 2 ^ + V(x)^-n-L4>, V(0)=Vo(z), 

where x G M d , for d = 2 or d = 3, respectively, and a < i.e. the nonlinearity 
is assumed to be energy- subcritical. In d = 2 the rotation term simply reads 

(2.2) fl-L=-i\Q\(x 1 d X2 -x 2 d Xl ). 

Remark 2.1. In d — 3 we could, without restriction of generality, choose a refer- 
ence frame such that fl = (0,0, |^|) T , yielding the same formula as in (2.2). For 
the sake of generality we shall not do so but consider the term il-L = — ffl'(iAV) 
in full generality. 

A potential V(x) G E d , satisfying (fi ■ L)V{x) =0, Va; G R d , is said to be 
axially symmetric (with respect to the rotation axis fl G M 3 ). In particular, this 
holds in the case of an isotropic trap potential, i.e. a potential of the form (1.2) 
with 71 = 72 = 73. 



CAUCHY PROBLEM FOR NLS WITH ROTATION TERM 



3 



Formally, (2.1) preserves the total mass 

M:= f \ip{t,x)\ 2 dx, 

and the energy 

(2.3) En~ [ hvyj\ 2 + V(xM 2 + ^—\yj\ 2 °+ 2 -^(n-L)4>dx. 

Note that, the last term is indeed real valued (as can be seen by a partial inte- 
gration). In order for these two quantities to be well defined, we shall study the 
Cauchy problem corresponding to (2.1) in the space 

E := {/ G if 1 ^) : \x\f G L 2 (R d )}, 

equipped with the norm 

11/11=:= ll/ll| 2 + ||V/||| 2 + |kf||| 2 . 

We remark that even if the potential V(x) is chosen to be identically zero, it would 
not be enough to consider the Cauchy problem for tj) G H 1 (R d ), since in this case 
we can no longer guarantee that 

(2.4) L a := ip(Sl ■ L)%l)dx < +oo. 

Physically speaking, this means that ip has finite angular momentum. The choice 
of E is therefore natural in our situation and not necessarily linked to the presence 
of a harmonic trapping potential, in contrast to [5, 6, 7]. 
We can now state the main result of this work. 

Theorem 2.2. Let < a < 2/(d - 2), A e 1, 11 e R d , for d = 2, 3 and denote the 
smallest trap frequency by 7 := min{7j}^ =1 . 

(1) Then, for any given initial data ipo G E, there exists a unique global in-time 
solution tp G C([0, 00); E) to (2.1), provided one of the following conditions is sat- 
isfied: 

(i) the nonlinearity is L 2 -subcritical a <2/d, or 

(ii) a ^ 2/d and A ^ 0, i.e. the nonlinearity is defocusing. 

(2) On the other hand, if A < 0, and if either: 

(i) (Cl ■ L)V = 0, i.e. V is axially symmetric, and a ^ 2/d, or 

(ii) (Q • L)V 7^ 0, |Q| ^ 7, and a > a^/d, where 

(2.5) a n := 

Then there exist initial data ipo G E such that finite time blow-up of the correspond- 
ing solution ip(t) occurs. 

In fact, we shall prove Assertions (l)(i) and (ii) under the more general assump- 
tions on V(x), see Assumption 3.1 below. This, together with the fact that no 
condition on the size of |0| or {jj} d =1 is required, generalizes the earlier results 
given in [11, 12]. 

Concerning the possibility of finite time blow-up, we see that one has to dis- 
tinguish between the case of axially symmetric potential and the case where this 
symmetry is broken. The reason will become clear in the proof given below. In 
the case of a non-axially symmetric potential we can rigorously prove the occur- 
rence of blow-up only in under the additional restrictions |f2| ^ 7, and a an/d, 
i.e. only for a limited range of nonlinearities. It is easily seen that in d — 3, the 
set of cr's satisfying our conditions is non-empty, provided |f2| 2 < I7 2 . Also note 
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that in the case of vanishing rotation lim|^| — a n = 2, yielding the usual range 
of L 2 -supercritical nonlinearities. At this point it is not clear if these additional 
restrictions are only due to the strategy of our proof, or if they indicate an actual 
difference in the behavior of solutions to (2.1). In particular, the question whether 
or not finite time blow-up occurs in situations where SI > -7 is completely open so 
far. In terms of physics, the latter would correspond to the case where the rotation 
is stronger than the trap and thus one would expect a behavior which is similar (at 
least qualitatively) to the "free" case, i.e. without any potential. 

This paper is now organized as follows: Section 3 is devoted to the proof of 
Assertion (1) of Theorem 2.2. To this end, we shall first prove local in-time existence 
for solutions to (2.1). Also, we shall see that a naive use of the conservation laws 
for mass and energy in general leads to restrictions on |f2| or {■Jj}j—i- We shall 
show in a second step how to overcome this problems using a coordinate-change. 
Assertion (2) of our main Theorem is then proved in Section 4 and we finally collect 
some concluding remarks in Section 5. 

3. Local and global existence 

In this section we shall allow for more general class of potentials V(x) satisfying 
the following assumption. 

Assumption 3.1. The potential V : K d — s> R is assumed to be smooth and sub- 
quadratic, i.e. for all multi-indices k E N d , with |fc| > 2, there exists a constant 
C = C(k) > such that 

(3.1) \d k V(x)\^C forallxeM d . 

Remark 3.2. Clearly, a harmonic trapping potential of the form (1.2) is sub- 
quadratic. Assumption 3.1 allows us to take into account more general situations 
of physical interest, such as a combined harmonic trap plus optical lattice potential, 
see e.g. [9]. Note however, that under Assumption 3.1, the potential is not neces- 
sarily bounded below (or confining). In particular we can also allow for repulsive 
potentials such as V(x) = — 7 2 |x| 2 , see [5]. 

As a first, preliminary step, we shall prove the following local well-posedness 
result. 

Lemma 3.3. Let ipo e S, w 6 R, and < a < 2/(d — 2). Moreover, assume that 
V satisfies Assumption 3.1. Then there exists a time T = T(||^o||s) > and a 
unique maximal solution ij) <E C([0,T);E) of equation (1.1) with tp(0) = ipo- The 
solution is maximal in the sense that, if T < +00, then 

lim \\Vip{t)\\ L 2 = +00. 

t^rT 

Moreover, the following conservation laws hold 

(3.2) M(i) = M(0), Ea(t) = En(0), 
whereas for the angular momentum we have 

(3.3) Ln(t)+ [ f i\ip\ 2 (fl ■ L)V(x)dx = Lq(0). 

JO JR d 

Proof. The proof is an adaptation of classical arguments, based on a contraction 
mapping (via Duhamel's formula) and Strichartz estimates for the linear (unitary) 
group U(t) — e' ltH generated by the Hamiltonian 

H = ~A + V(x)-n-L. 
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In our case, Strichartz estimates can be obtained by following the approach of 
Kitada [14]. All we need to do is to check the assumptions given there: First, we 
note that the Hamiltonian is a quantization of the classical phase space function 
H(x,t;) = — 1£| 2 /2 + V(x) + id ■ (x A £), which is smooth and sub-quadratic in x 
and £, due to Assumption 3.1. We easily check that also the second assumption of 
[14] holds true, i.e. 

{H{x^),e^W)m) = (H(x, V ),e- ix -\(x)Jtn)), 

where (•, •) denotes the usual duality bracket between tempered distributions and 
Schwartz functions. This expression is equivalent to the fact that H is essentially 
self-adjoint. We therefore conclude that there is a 5 > such that 

In particular it follows that the Strichartz estimates for H are analogous to those 
found in the well-known case of NLS with quadratic potentials [6], i.e. the rotation 
term does not influence the dispersive behavior (locally in time). The existence of 
a local in-time solution then follows analogously as given therein. The conservation 
laws (3.2) follow from straightforward calculations in combination with a standard 
density argument (see e.g. [8]). Finally, in order to prove the blow-up alternative 
we first compute 

(3.4) ^-\\xm\\h = 21m / xm^m & < w^mwh + \\vm\\l>- 

Thus, as long as || VV'(i) ||z, 2 is bounded, Gronwall's inequality yields a bound on 
||a;^(t)||£2 as well. In view of mass conservation, the only obstruction to global 
existence is therefore given by the possible unboundedness of ||VV>(*)||l 2 m [O,? 1 ]. 

□ 

Remark 3.4. For quadratic potentials of the form (1.2), Strichartz estimates can 
be obtained explicitly by invoking a generalization of Mehler's formula for the kernel 
of U(t), c.f. [7]. Indeed, by making the following ansatz 

where Hj{t) 6 R+ and F(t,x,y) is a general quadratic form in x and y with (yet 
to be determined) time-dependent coefficients. Substituting this into the linear 
Schrodinger equation yields a coupled system of differential equations for these co- 
efficients. Solving this system, however, is in general rather tedious. This approach 
is therefore only feasible under some simplifying assumptions, such as Q — [5, 6], 
or V(x) = ^-|x| 2 with |Q| = 7 as it is done in [11, 12]. 

In view, of (3.3), we immediately conclude the following important corollary. 

Corollary 3.5. IfV(x) is such that (O • L)V — 0, then we also have conservation 
of angular momentum, i.e. Ln(t) = Lq{0), and in addition it holds 

(3.5) Eo{t)= [ hvyj\ 2 + V(xM 2 + ^-\yj\ 2 ° +2 dx = E (0). 

Thus, in the case of axially symmetric potentials V(x), there are in fact two 
conserved energy functionals corresponding to (2.1). 

With a local existence result in hand, we can ask about global existence. In order 
to infer T = +00, one usually invokes the conservation of mass and energy (3.2). 
The problem is, that due to the appearance of the angular momentum rotation 
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term, the energy En(t) has no definite sign even if V ^ and A ^ (defocusing 
nonlinearity) . A possible strategy to overcome this problem is to rewrite the linear 
Hamiltonian as 

1 1 I o 1 2 

(3.6) H = --A-Q-L + V(x) = -(-iV - A(x)f + V(x) - Lj-r 2 , 

where A(x) — ft A x and r = a; A ft|/|ft denotes the radial distance perpenticular 
to ft. Note that A(x) can be considered as the vector potential corresponding to a 
constant magnetic field B = VA A = 20. The corresponding "magnetic derivative" 
Da '■— — i(V + A(x)) is known to satisfy, cf. [8, Chapter 7]: 

It can therefore be used to control the nonlinear potential energy oc ||^|| L 2<t+2 via 
Gagliardo-Nirenberg type inequalities. If in addition, V{x) is given by (1.2) with 

|ft| ^ 7 we infer that V(x) — ^-r 2 > 0. In this case, the linear part of the energy 
is seen to be a sum of non-negative terms, and global existence can be concluded 
as in the case of NLS with quadratic confinement [6]. However, it seems impossible 
to extend this approach to situations in which |ft| > 7, even if A > 0. In order to 
do so, we shall follow a different idea, which invokes a time-dependent change of 
coordinates. 

Proof of Assertion (1) of Theorem 2.2. We start with the i 2 -subcritical case, i.e. 
< a < 2/d which follows by standard arguments. Namely, we write the nonlinear 
Schrodinger equation as the solution of a fixed point equation, using Duhamel's 
formula: 

1>(t) = U(t)il>o-i\ J U{t- S )(\^\ 2a { s )^(s)) ds 

=■■ *W0(t). 

Next, we calculate the commutators [V, H] and [x,H]. Explicitly, we find 

[V, H] = - 1[V, A] + [V, V] + i[V, SI ■ (x A V)] 

= \7V + i[V,x- (V Aft)] 
= W + iV AQ 

by the well-known formula a ■ (b A c) = det(a, b, c) = (a A b) ■ c for three-dimensional 
vectors a, 6, c. Similar calculations yield 

[x, H] = V - ifl A x 

From this we can deduce 

V$(V>)(*) = U{t)Vip - iX [ U(t - t)V (|^(r)| 2ff V(T)) dr 

Jo 

-if U(t-r) (VV - ifl A V) *(V>)(r) dr, 
Jo 

and 

x$(ip){t) = U(t)Vip -i\ f U(t-r)V (\xP(T)\ 2a iP{T)) dr 

Jo 

-if U(t - t) (V - ifl A x) $(^)(r) dr. 
Jo 

Using Strichartz estimates (see the discussion in the proof of Lemma 3.3) we have 
that 

IIV ; lkp(0,T;L9)nL~(0 1 T;L 2 ) < Cll^olU 2 + C 1 1 "0 1 1 Z,fe (0,T; i9 ) , 
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where 



4cr + 4 2cr(2cr + 2) 

^ 2ff + 2 ' fc= 2-(d-2) g - 



If a < 2/d it holds that 1/p < l/k and thus 

IMIl*(0,T;£«) ^T^-^UW^T.^y 

If we choose T — T* > small enough, we can absorb the last term on the right 
hand side in order to get a bound on ||V'||z>(o,T*;L<!)nL°°(o,T*;L 2 )- Since we can 
shift the time interval [0,T*] by an arbitrary amount of time, in the same way 
we can get a uniform bound on \\ip\\Lp(i:Li)nL°°(i:L 2 ) f° r every interval of length 
|7| ^ T*. Thus, by splitting any arbitrarily large time interval [0,T] into suffi- 
ciently small sub-intervals {I n }n=i sucn that \I n \ ^ T* and iterating the bound 
\\ip\\LP(i;L<i)nL°°(i;L2) < C* , wc i nfcr II ip 1 1 Lp(o,T ; L<i)nL~ (o,t ; i, 2 ) < C where C < +oo 
depends on the value of T. Moreover, since ip £ E we also get 

11^11^(0, T;L-!)nL~(0 : T;L 2 ) + l|VV'|Up(o,T;L'!)nL~(0 : T;L 2 ) 
< CUVolls + C||V||i < L (0l T;L«) (IWIIl»-(0,T ; L«) + ||VV||LP(0,TiL«)) + 

+ CT (||o;V|Uoc ( o,t;l 2 ) + II V^|| l ~ (0 ,t ; l 2 )) 

From here, we proceed as before to obtain a uniform bound for the left hand side 
for small T* and thus by iteration for arbitrary time intervals [0,T]. 

Next we consider the case of an L 2 -supercritical nonlinearity a > 2/d. In this 
case, the iterative argument given above breaks down. The basic idea is to use a 
change of coordinates in order to bring equation (2.1) into a more suitable form. 
For the sake of notation we shall only consider the case d = 3 in the following. We 
first note that by using the skew- symmetric matrix 

/ fl 3 -ft 2 N 
6 := -fi 3 fti 
\ tt 2 -Oi 

the wedge product with the angular momentum can be written as 

A x = -9 • x. 

Then the matrix-exponential 

X(t,x) :=e et -x 

defines a rotation of the vector x £ M? around the axis O by an angle of — jO|i. Its 
time-derivative can be calculated as 

(3.7) d t X(t, x) = e- X(t, x) = —Q A X(t, x). 

Denoting the wave function in rotated coordinates via 

$(t,x)=ip(t,X(t,x)), 

we conclude from (3.7) that 

id t i){t, x) = id t ip(t, X(t, x)) - i (Q. A X(t, x)) ■ Vip(t, X(t, x)). 
Rewriting — i(Q A X) ■ V = — iQ ■ (X A V) = Cl ■ L, we arrive at 

idt$ = + am 2 > + v(x{t, x ))$. 

where we have also used the fact that the Laplace operator is invariant with respect 
to rotations, i.e. 

A x1 p(t,X(t,x)) = A x ^(t,X(t,x)). 
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Dropping all the tildes and denoting W(t, x) = V(X(t,x)), we conclude that up 
to a change of coordinates, equation (2.1) is equivalent to the following NLS with 
time-dependent potential 

(3.8) id t ip = -^Aip + X\ip\ 2a ip + W{t,x)i) 

Note that W(t,x) is smooth w.r.t. iel and sub-quadratic w.r.t. x 6 M 3 with the 
same (uniform) constants C{k) as given in Assumption 3.1 for V(x). Moreover, if 
V{x) is axially symmetric, i.e. (il ■ L)V(x) = 0, equation (3.7) implies that 

d t W(t, x) = -ft A X(t, x) ■ VV(X(t, x)) = -i(fi • L)V(X(t, x)) = 0, 

and hence Wit, x) = W(0, x) = V(x). The energy corresponding to the transformed 
NLS (3.8) is given by 

Ew{t) := J ^\ViP(t,x)\ 2 + \m,x)\ 2 °+ 2 + W(t,x)m,x)\ 2 dx. 

However, since the potential W{t,x) in general is time-dependent, the Ey/{t) is no 
longer a conserved quantity. Rather, we obtain that 

(3.9) J t Ew(t) = j W^MM)) 2 dx - 

Nevertheless it is not hard to prove Assertion (l)(h) of Theorem 2.2: In view of the 
blow-up alternative, stated in Lemma 3.3, it suffices to show || Vip(t)\\ L 2 < +oo, for 
all T > 0. To this end, we first estimate 



±\W(t)\\h <E W ® + 



J W{t,x)\iP(t,x)\ 2 dx 



^E w (t) + C\\x^(t)\\ 2 L2 , 



under the assumption that A > 0. Integrating equation (3.9) and having in mind 
that W(t, x) is sub-quadratic in x, we obtain that 

(3.10) \\Vm\\l'<E w {0) + J ^-E w ( s )d s + C\\x^(t)\\h 

<c (i+\\x^(t)\\b+j* iN>(s)in 2 d s y 

Recalling inequality (3.4), we infer 



d_ 
~dt 



Mt)||i a + < C (l + \\xm\\h + f \\^{s)\\h ds 



which by Gronwall's inequality yields an uniform bound on ||xV(£) \\l 2 for every 
time interval [0,T]. With this in hand, we can bound ||VV>(i)||L 2 by simply using 
inequality (3.10) once more. □ 

Remark 3.6. In particular, for ft = (0,0, |ft|) T and V(x) given by (1.2), we 
explicitly find 

W(t, x) = i( ( 7l 2 cos 2 (|0|i) + 7 2 sin 2 (|ft|i)) x\ + ( 7l 2 sin 2 (|ft|t) + 7 2 cos 2 (|ft|t)) x\ 

+ sin(2|ft|i) ( 7 i - ll) x lX2 + j 2 xfj . 
Clearly, W — 5(71 a: 2 + ~j 2 x\ + 73X3) in the axially symmetric case 7 2 = 7! • 
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4. Finite time blow-up 

This section is devoted to the proof of assertion (2) of Theorem 2.2. It follows 
from the following lemma. 

Lemma 4.1. Let A < 0, a < 2/(d -2), fie M d , for d = 2,3, and V(x) be a 
quadratic potential of the form (1.2). Denote 7 = min{7j}^ =1 and let an be as in 
(2.5). If either 

(i) (ft • L)V = 0, cr ^ 2/d, and £ (0) < 0, or 

(ii) (ft • L)V + 0, |ft| < 7, cr ^ an/d, and Eh(0) < 0, 

iften </ie corresponding solution necessarily blows up in finite time. 

Note that the condition for the energy of the initial data tpo are not identical in 
both cases. The reason will become clear in the proof given below. 

Proof. To simplify the arguments later on, let us first compute the conservation 
laws for the mass and momentum densities, i.e. p := \ij)\ 2 and J := Im(r/>V^>). 
Indeed a straightforward calculation yields 

(4.1) dtp + div J = iil ■ hp. 

On the other hand, for the current density J we find 



(4.2) 



d t (im^W)) = Im ^ (-^AV> + iV(x)^ + iX\^\ 2 iP + in ■ Lip^j W 

+ Im (Vv (^Atp - iV(x)ip - iA|Vr _1 V + ■ ) 

Next, we calculate 

Im (^V(ift • LV)) = Im (^(ift • i)V^) — ft A J. 

Thus we can combine the two terms in (4.2) which stem from the rotation via 

Im ((in ■ L)lpVil)) + Im (lp~(in ■ L)Vtp) - ft A J = (in ■ L) J - ft A J 

where we have used that ift • L is real-valued. The other terms in (4.2) are usual 
in quantum hydrodynamics, see e.g. [2], yielding the following equation for J: 

(4.3) d t J+div(Re(V^O Vii))+-^V\ip\ 2a+2 +pVV = ^AVp+(ift-i) J-ftA J. 

The proof of finite time bow-up now follows by the classical argument of Glassey 
(see [8]). To this end, we consider the time evolution of 

m--= \ [ \x\ 2 m, X )\ 2 dx. 

z JR d 

Differentiating with respect to time and using (4.1), we obtain 

^-I(t)= f x-J(t,x)dx+ f ^-(ift • L)p(t,x)dx. 
dt J R d J R d 2 

Integrating by parts and using (ft • L)\x\ 2 — shows that the second integral in fact 
vanishes, i.e. we have 

d 



f 

-I(t) = x- J dx. 



dt 

Differentiating in time once more and using (4.3), we obtain 

j t J x ■ J dx = J x ■ ^ - div (Rc(V^ ® VV>)) - A— ^-y V|t/>| 2<t+2 - pVV + J AVp 
+ (in ■ L)J - n A J)dx, 
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which we rewrite as 

d f f d 

— / x- Jdx= / \V%1)\ 2 + A— ^—\^\ 2a+2 - px ■ W + x ■ (iSl ■ L)J 
dt J J a + 1 



(4.4) dt . 



x 



■Q A J dx. 



Now we first note that for any potential V(x) of the form (1.2) we have x- W = 2V. 
Moreover, we compute 



(O • (x A V)x) • J dx 

R d 



/ x ■ (id ■ LJ) dx = — I (fi ■ Lx) ■ J dx = — I 

JR d JR d JR 

= -/ ((fi Ai) • V)x ■ J dx = - / (CI Ax) ■ J dx, 

JR d JR d 



shows that the last two terms in (4.4) cancel each other. In summary we arrive at 
the following identity 

(4.5) ^I(t)= J m 2 + \^- 1 \^+ 2 -2V\^dx, 

which is in fact exactly the same as in the case of NLS without rotation, c.f. [8]. 

We can now prove assertion (i): Recall from Corollary 3.5 that if the potential 
V(x) is axially symmetric, then E (t) = E (0), with E defined in (3.5). Hence 
from (4.5) and V > we can write 

P I(t) 2E + X^—^ f \^\ 2a+2 dx 



dt 2 v ' " u a + 1 . _ 
Assuming E < 0, A < 0, and a ^ 2/d, we consequently obtain 

^ < ~ C > 

for some constant C > 0. Integrating this relation twice, we obtain 

C 

I(t) < --t 2 + Cl t + c 2 

with some integration constants c\ and C2- Thus, if the solution ip(t) G £ were to 
exist for all times, there would be a time T* < +oo, such that I(T*) < 0. This 
however is in contradiction with the fact that, by definition, I(t) ^ for all t € K 
and hence the assertion is proved. 

In order to prove assertion (ii) we again consider (4.5): The problem is that in 
the case of a non-axially symmetric potential (fl ■ LV(x) ^ 0), the energy E is no 
longer conserved. Rather we only have the conservation law for Ea(t) = Eq(0). In 
order to use this piece of information, we first add and subtract to (4.5) a multiple 
of the angular momentum La(t), i.e. 

4*i{t)= I |v?M 2 + ^-|V'l 2CT+2 -2y|v| 2 -a^-i^dx+ [alffl-L^dx, 

dt 2 J Rd a + 1 J 

where a > is a parameter to be chosen later on. Using Cauchy-Schwarz and 
Young's inequality the last term on the r.h.s. can be bounded by 



JR 



R d 1 M 



where 6 > is another free parameter to be chosen later on. We consequently 
estimate 



/ aipQ-Ltpdx. 

JR d 
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Now, we choose such that 2(1 + ^) = a, that is 9 = 9L ^-- In this way we have 

+ [ \^-^m 2 ° +2 dx+ [ (-( a + 2)V+^-\xA\^dx. 
Let 7 := min(7i, 72, 73), and choose a such that 



(a + 2) 7 2 = a2 l^l 2 



2(a - 2)' 
This yields a = an with 



an 



4 7 2 



By doing so, the last term in the previous inequality is seen to be non-positive and 
furthermore we conclude that, for A < and a > ^p, it holds: 

d 2 

(4.6) -pj(t) < an^n(t) = a n E n (0). 

Thus, if the initial energy Eq(0) < the second derivative of I(t) is again negative 
and we can argue (by contradiction) as before. □ 



5. Concluding remarks 

As we have seen above, equation (2.1) can be considered (upon a change of coor- 
dinates) as a special case of NLS with time-dependent potentials (sub-quadratic in 
x) . This class of models has recently been studied in [7] . Following the arguments 
given therein, one could infer global in-time existence of (2.1) for sufficiently small 
initial data -00 £ regardless of the sign of the nonlinearity. Moreover, growth 
rates for higher order (weighted) Sobolev norms can also be obtained as in [7]. In 
addition, we note that for a repulsive, isotropic quadratic potential V(x) = — ^\x\ 2 , 
the timc-dependent change of coordinates is trivial and we could henceforth con- 
clude global in-time existence for sufficiently large 7 > by following the arguments 
given in [5]. 

We also want to point out that for the usual NLS with u = 2/d there is an extra 
symmetry which has been successfully deployed in the study of blow-up (yielding 
explicit blow-up solutions and blow-up rates), see e.g. [22]. Using the so-called 
Lens transform [13] one can transfer (most of) these results to the case of NLS with 
isotropic time dependent quadratic potential W(t,x) = 7(t)|x| 2 , see [7]. However, 
it is argued in [7] that such an approach is only feasible in the case of isotropic 
potentials and thus, we cannot expect from it any further insight on the possibility 
of blow-up in our case, when (L ■ Q)V(x) ^ and |f2| > 7. 

Finally, it is worth noting that the effect of the angular momentum rotation 
term in our model is very different from other situations. For example, it has 
been shown for the Euler equations with Coriolis force that blow-up can be delayed 
through a sufficiently strong rotation term [18] (see also as [3] for a related result). 
Clearly, the situation in our model is much more involved, and we can not expect 
an analogous result to be true (the counterexample being the case where V(x) is 
axially symmetric). 
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